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The aviation sector is profitable, but sensitive to economic fluctuations, geopolitical constraints and govern- 
mental regulations. As for other means of transportation, the relation between origin and destination results in 
a complex map of routes, which can be complemented by information associated to the routes themselves, for 
instance, frequency, traffic load or distance. The theory of networks provides a natural framework to investigate 
dynamics on the resulting structure. Here, we investigate the structure and evolution of the Brazilian Airport 
i£>s" ' Network (BAN) for several quantities: routes, connections, passengers and cargo. Some structural features are 

j.-^ , in accordance with previous results of other airport networks. The analysis of the evolution of the BAN shows 

*vj ' that its structure is dynamic, with changes in the relative relevance of some airports and routes. The results 

indicate that the connections converge to specific routes. The network shrinks at the route level but grows in 
Sv , number of passengers and amount of cargo, which more than doubled during the period studied. 
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Ph! I. INTRODUCTION 

6 . 

O , During the 20th century, society witnessed a large change in terms of transportation options. In the first half of the century, 

. transportation was dominated by ships in intercontinental routes while trains were the most important carriers in continental 

^ routes. The aviation sector developed at a fast pace during these years, nowadays it dominates intercontinental traffic and it has 

. ,— I a strategic importance in continental flows. This sector has an increasing participation in the domestic transportation market, 

^^1 competing with other more traditional transportation categories. The short travel time is the major motivation to choose air 

^-T *. transportation, even though the flight prices are still high in some regions. The flight velocity together with long range routes 

O have deep implications in terms of transportation. People can move fast between different regions and consequently, spread 

I— I information and diseases faster as well. It was already observed that the world wide airport network is of importance for global 

epidemics j^. 

The development of the air transportation naturally results in a complex map of routes established between airports. These 
routes can be regular or non-regular and they are mainly constrained by economic, geopolitical, regulatory and technological 
QQ \ factors. This cornplex structure of flights can be mapped into a network and analyzed by using tools from the so-called cornplex 
f^ ■ network theory Ij, B m^ Indeed, the airports network structure has been investigated before on the world wide HI 3, S la, 
C^ ' U] and domestic levels Q IH H U^ Hll ll2. U^ I3l^ involving countries with different economic and politic situations, and 
T-j- . population/area sizes. Some attempts to model the world airport network (WAN) have been proposed considering geopolitical 
f^ ' constraints 01511 . passengers behavior 01611 and optimization principles lUTIl . 

OO In particular, the WAN has skewed distributions for degree, passengers traffic and betweenness centrality in an extended 

^^ range |5, 6, 7]. In a weighted version, BaiTat and collaborators found a strong coiTelation between the number of routes and the 

►>. passenger traffic in an airport, and a Unear relation between the average passengers traffic and the betweenness ^ |al- Using 

another database, Guimera et al. suggested that there are some cities with a small number of routes but large betweenness which 

could be related to geopolitical constraints [7]. Moreover, they found a hierarchical organization where less connected nodes 

^ generally belong to well interconnected communities, while hubs connect many airports that are not directly connected 1 7 ] . In the 

■ - ' case of the Chinese and US airport networks, different attempts were made to fit the degree distribution which in principle, obeys 

double Pareto's law |la,|M,[i2l- On the other side, the Indian and Austrian networks apparently follow a single power-law Ill0l[l4ll . 

Although time evolution analysis has been already performed in a short time-scale basis [jSl |3, u^ llM and shows that there 

are only small fluctuations in some structural features of the airport network during the week, the time evolution analysis in a 

year-scale has not been performed yet. The importance of such study is related to the long term evolution of the airport network 

structure. Such structural evolution reflects the economic health of the system (air transportation related companies and the 

country itself), since routes are inserted and removed according to passengers and cargo demand, and new airports are built and 

expanded according to regional development. The network evolution is also constrained by governmental policies which regulate 

the sector 1 18] and by economic fluctuations (e.g. oil prices) which directly affect transportation sectors such as aviation i 1 8|] . 
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The following analysis of the evolution of the Brazilian airport network uncovered a number of features, suggesting that de- 
spite of small fluctuations in some structural properties, the network has undergone a strong rewiring dynamics in the considered 
period. What in principle might appear as a slowdown of the aviation sector, in fact is a reorganization of the structure according 
to changing demand. One of the consequences has been an increase in the centrality of some airports, which might affect the 
performance of the network. 

II. THE BRAZILIAN AIRPORT NETWORK 

Extensive literature on the fundamentals of complex networks is available IS S Hi. however, some basic definitions will 
be shortly introduced before the discussion of the results. A network F := {V,E) can be defined as a set V of nodes i 
{i — 1,2, ...,N) and a set E of K edges {i, j) connecting these nodes according to some rule. Another set of real numbers 
W is also defined, such that one or more elements of this set can be mapped into each edge as weights, in order to represent a 
specific feature of the system. Representing the network as a matrix A (adjacency matrix), the element a^ is 1 if there is an edge 
connecting nodes i and j (the direction of the edges can be distinguishable). Another matrix W (weight matrix) can be defined 
such that each element carries information about the weights assigned to the respective edges. By these definitions, the degree 
(total number of edges or direct neighbors) and the strength (total weight of the edges) of anode i are respectively k{i) = ^ a^ 
and s{i) — ^ aijWtj . In a directed network, not all edges are bidirectional, which means that {i,j) is not necessarily the same 
as (j, i). To quantify this asymmetry, we define the reciprocity R to count the proportion of bidirectional edges in the network 
(eqn.[T}. 

R= ^^\ -± (1) 

A path of size d connecting nodes i and j is defined as a sequence of d — 1 nodes from the set V such that there is an edge 
connecting the node m with the node m + 1 of the sequence, being the first and last nodes in the sequence, respectively, i and 
j. If there is a path between % and j, it means that j is reachable from i. A geodesic of size d {af^ is the shortest path between 
nodes i and j. The diameter Ddir is defined as the largest geodesic in the network, which is the minimum distance necessary to 
reach any pair of nodes. The average distance ddir measures the average minimum distance between any reachable pair of nodes 
(eqn.|2]i. 

(^-)-M]73T)E4 (2) 

We can measure the amount of geodesies going through a specific node by using the Freeman betweenness centrality \\y^- 
Given that | cr,;j (K) \ is the number of geodesies between nodes i and i that pass through node h and | Uij \ is the total number 
of geodesies between i and j, we get the proportion of shortest paths going through specific nodes by using equation|3] 

The Brazilian airport network (BAN) is so defined by mapping Brazilian cities (henceforth referred only as airports unless 
specified) into nodes i and associating a directed edge [i,]) between two nodes i and j, whether the respective cities are served 
by airports containing non-stop regular flights (henceforth named route [i, j)) (see fig. [T]). By this definition all the non-regular 
flights (e.g. air charters, private, military) as well as airports without regular flights are discarded. For each year, from 1995 to 
2006, a multi-layered network BAN year with Nyear nodes and Ky^ar edges is constructed. Each layer has the same set of nodes 
and edges, but represents different features associated to the respective route and represents different weights luOll . In the present 
case, the network is defined with 3 layers: the first (L^) is associated to the number of connections (or flights) completed during 
the year in the respective route. The second (L^) corresponds to the number of passengers and the third (L^) to the amount of 
cargo (in kilograms) carried by these flights. In this specific system, the multi-layered approach is equivalent to assigning three 
different weights for each edge in the network. 
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Figure 1: (Color online) Routes starting/ending in Sao Paulo city in 2006. Some important cities are highlighted, e.g. the capital Brasilia. 
Support material includes the visual evolution of these routes from/to Sao Paulo along the years. 

The data used to build the BAN are publicly available in the web-page of the Agenda Brasileira de Aviagao Civil (Brazil- 
ian National Agency of Civil Aviation) and they are directly extracted from the electronic files of the Anudrio Estatistico do 
Transporte Aereo (Air Transportation Annual Report) which contains specific information about flights and statistics about the 
domestic and non-domestic airline companies operating in Brazil |21]. The airports are identified either by the name of the 
main served city ( 75%) or by the airport name (^ 20%). Approximately 5% of the airports serve the same metropolitan area, 
i.e. are settled within the same city or nearby city. Airports within the same metropolitan area are considered as one node in 
the present network. The data are filtered to remove self-loops, which in principle are meaningless and appear specially in old 
files for unknown reasons (possibly misprints). Multiple entries are also removed and the respective data are summed into a 
single entry (except in the cases where misprints are clearly identified). It is important to note that flight information of regional 
and national companies are separated till 1999, these data are also summed in order to represent the complete national flight 
information. All the networks are observed to be completely connected with the exception of one year (1999), which besides the 
giant component, contains 2 small connected components with 2 nodes each. This feature has negligible effects in the following 
analysis and we can safely assume that all networks are connected. 

For comparative analysis, the network is randomized maintaining the degree distribution \^ with self-loops and multiple 
edges being forbidden. Observe that bidirectional edges are not considered as multiple edges. The measurements are obtained 
by averaging over 10 network samples. This randomization mechanism is adopted to skip the influence of the degree distribution 
in the results. 

III. STRUCTURE AND EVOLUTION 
A. Unweighted BAN 

The configuration of the BAN is observed to be dynamic. Despite some fluctuations, the total number of cities Nyear (i-e. 
nodes) and average number of routes ki„+out (i-C- directed edges) decreased consistently as observed in Table U For the sake 
of comparison with previous results of other airport networks, the average degree (fc„„d) and clustering coefficient (ccund) Hi 
of the undirected version of the BAN are also presented in Table |T] The clustering coefficient measures the probability of two 
common neighbors of a node i being connected with each other While the average degree decreased by approximately 28.3% 
in 2006 if compared with 1995, the values of the clustering coefficient are relatively high and stable during the period. In 
comparison to other airport networks, the BAN is positioned as the second largest national network after the USA when ranked 
by the number of airports J^, uM- Moreover, the high clustering coefficient and the small average shortest path length (ddir) 
(eqn.[Tli indicate that the BAN is small-world, which is in accordance with the other national networks. 

In previous papers, the authors usually suggest that with good approximation the direction of the routes are irrelevant 151 la 
0>H[l3l for those networks. Nevertheless, by measuring the reciprocity R (eqns.[T]l of the BAN, it is clear that although highly 
reciprocal, R is not close to 1 (Table|I]l, which means that one can fly from one airport to another but not necessarily return using 
the same route. Therefore, to maintain this path feature of the network, hereafter the directed version of the network will be 
adopted. 



Year 


1995 


1996 


1997 


1998 


1999 


2000 


2001 


2002 


2003 


2004 


2005 


2006 


N 


211 


224 


234 


219 


208 


183 


171 


172 


143 


152 


154 


142 


{f^und) 


13.19 


11.49 


10.63 


10.49 


10.46 


10.60 


11.09 


9.79 


9.71 


9.49 


9.73 


10.28 


{cCund) 


0.66 


0.65 


0.62 


0.63 


0.66 


0.66 


0.62 


0.63 


0.62 


0.61 


0.65 


0.63 


\r^in-\-out) 


22.84 


20.37 


18.81 


18.80 


18.86 


19.04 


19.64 


17.70 


17.55 


17.16 


17.65 


18.32 


R 


0.836 


0.865 


0.864 


0.879 


0.885 


0.880 


0.863 


0.888 


0.887 


0.888 


0.891 


0.869 


Dd„ 


6 


5 


6 


6 


6 


5 


5 


6 


6 


5 


5 


5 


{ddir) 


2.44 


2.47 


2.36 


2.39 


2.34 


2.49 


2.40 


2.39 


2.36 


2.26 


2.36 


2.34 


Jin 


3.10 


3.09 


3.05 


3.04 


3.01 


2.91 


3.00 


2.93 


2.92 


2.87 


2.91 


2.92 



Table I: Evolution of the basic structural features of the BAN network. 



The average shortest path length (ddir) (eqn- EJ is in general, slightly higher ('^ 5%) than the randomized version of the 
same network, while the diameter D^ir is about 12% larger These results might be related to geographical constraints which 
are unavoidable in the BAN due to the size and the shape of the country, some small (regional) airports are only connected to 
closer airports and have no long-range routes due to restricted traffic demand. The same effect could explains why the diameter 
fluctuated little during the period (Table |I|. 

The high correlation (Pearson coefficient r ^ 0.99) between the number of in- and out-degree permit us to analyze the 
in-degree distribution and extend the results to the out-degree. Since the BAN is relatively small for statistical purposes, the 
cumulative data are taken to diminish the data spread, but even so, it is not conclusive which is the right functional form to best 
describe the in-degree cumulative distribution. A reasonable fit is obtained by using a stretched exponential P{k*in ^ kin) = 
C exp{—/3k~^) (Fig. 121). We verify that in recent years, some of the most connected airports are operating close or even over 
the estimated capacity, which means that possibly the decay is a result of physical constraints in the airports as pointed out 
before [5J. 
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Figure 2: (Color online) Cumulative in-degree distribution of the BAN for (a) 1995 and (b) 2006. Reasonable good fits are given by stretched 
exponential functions with parameters (a) C = 1.414, f3 = 0.3196 and a = 0.631, and (b) C = 1.648, /3 = 0.474 and a = 0.559. 

A better picture of the evolution of the degree heterogeneity (e.g. in-degree) is obtained by measuring the entropy Sdegree of 
the distribution Pdegree (fc) 12211 (eqn. |4]l. In this context, small entropy indicates that the network has a well known degree while 
high entropy show the opposite, i.e. a diversity in the degree values. The results (Table IJi suggest that the nodes have became 
less heterogeneous (lower entropy) with time. As pointed out before, a more homogeneous network is less robust to random 
failures [22]. A comparative analysis during the years suggest that the continuous skewed feature of the degree distribution 
maintain the network vulnerable to hub fails (e.g. congestion), but although small, the increasing homogeneity might contribute 
to make the network more fragile against random failures as well. 
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On previous analyses, the authors argue that the world [7] and the Chinese f§\ airport networks have a relatively stationary 
structure within a period of a year, however, the measurements fluctuations observed so far suggest that a rewiring dynamics 



is taking place over the years in the BAN. Investigating the establishment of new routes, we observe that a rewiring process 
underlies the evolution of the BAN, though routes connecting new airports 1 24] to the BAN also play an important role (Fig. [3]). 
Nearly 18% of the total routes in a year are new routes established betwen old airports (case aOO UM ). The number of added 
routes between old and new airports is less than 8% (aON). On the other hand, about 20% of the routes from the previous year 
are deleted between old airports (dOO), while the number of routes added between old and removed airports is just about 6% 
(dOR). The number of routes added between new airports (aNN) and deleted between removed airports (dRR) are irrelevant, 
approximately 0.2% and 0.1% (Fig.O. This result reflects the natural tendency of airline companies to optimize ground staff 
and use the same facilities rather than investing heavily on routes to new airports. 

Important to observe that only 100 airports are identified in all years but 395 different airports are counted during the period. 
Most of the airports appear and disappear during the years, but some of them stay in the network for a while and play an important 
role before being removed. The first case contains airports that temporary serve other nearby airports under makeover or even, 
unsuccessful tentatives of new routes. It also happens that small airports are included in the network for a short period of time. 
The airports in the second case are still important in terms of other services (e.g. air charter and private flights), however, due to 
several reasons, regular services are moved to specific cities in some regions. 

30 



20 



10 



-10 



-30 




nzi aOO 
^ aON 
^ dOO 
nZl dOR 
^ aNN/dRF 



1 9961 9971 9981 99920002001 20022003200420052006 
Year 



Figure 3: (Color online) Percentage of added (a) (based in the current year) and deleted (d) (based in the previous year) routes related to new 
(N), removed (R) and old (O) airports. Case aOO: added routes between old airports; aON: added routes between old and new airports; dOO: 
deleted routes between old airports; dOR: deleted routes between old and removed airports; aNN: added routes between new airports; and 
dRR: deleted routes between removed airports. 



This rewiring dynamics has important consequences in defining the paths between airports not directly connected. Since the 
average path length is larger than 1 and the diameter can be at most 6 (Table |I]i, one needs to go through intermediate airports 
to reach farther destinations and most important, sometimes the same akcraft and staff have to travel through many different 
airports during the day. 

The importance of the intermediate airports can be quantified by measuring the betweenness centrality (eqn.[3]l. The asym- 
metry in the relevance of some airports is observed on the (cumulative) betweenness distribution, which is relatively broad for 
all considered years. Figure |4] shows that in 2006 for instance, the cumulative distribution is reasonably approximated by the 
functional form P{B* ^ B) = Cexp(— /3i3") indicating that some airports are the most central ones. It suggests that most 
of the long range flights go through few intermediate airports which represent a socioeconomic relevance for a specific region 
and/or for the country itself. It is observed that the average betweenness has increased along the years (see SM - Fig. 2). This 
result indicates that flights between non directly connected airports are going through specific intermediate airports, meaning 
that some long non-profitable routes are disappearing and being replaced by routes connecting more important airports. This 
situation has drawbacks because the network becomes more dependent on these central nodes. Any problem related to delays 
(e.g. bad weather or holidays) on key airports can propagate and disturb the functioning of the whole network. Moreover, since 
most of the connections go through them, they should be effective for control of, for instance, smuggling, drugs traffic or disease 
spread. 

A path between two non-directly connected airports is not necessarily through the shortest path. Actually, many times al- 
ternative routes are available and normally, companies offer different travel options. Moreover, some companies (or even the 
passengers themselves) have the strategy of choosing cycles rather than going back and forth using the same sequence of airports. 
In a network, a cycle of size m is defined as a closed path that goes from/to node i, passing through to — 1 intermediate nodes, 
but not repeating a visited node. The BAN is observed to have more cycles than its random version with a relative peak in 1998 
for small cycles and in 2000 for larger ones (see SM - Table 1). 
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Figure 4: (Color online) Cumulative betweeness distribution of the real and random BAN in 2006. The data points are reasonable fitted by a 
stretched exponential with parameters C = 1.921, j3 = 6.629 and a = 0.254. 

B. Weighted BAN 

The routes themselves do not explain all the dynamics of the flights, specially, because the routes have different importance in 
terms of connecting different airports. If we consider the strength sj of a node i by summing the total number of connections (or 
flights - layer 1) during the year and divide it by the total number of routes (or (in+out)-degree), we can obtain the cumulative 
distributions in Figure |5] The resulting distributions corresponding to 3 different degree intervals are well described by the 
functional form Cexp{—(3x'^). The exponential parameter is larger for lower degrees and the 3 intervals overlap on some 
regions. These results indicate a non-trivial pattern of physical/traffic constraints according to the airport size. Airports with 
many different routes (the so-called hubs) have a lot of flights per route (at least 1.5 flights per day) while most of the low-degree 
airports have much fewer flights per route. The former case correspond to the busiest airports in the BAN, while the last case 
includes small airports operating few regular flights (or routes not operated in a daily-basis but in some specific days during the 
week, which is due to several reasons, e.g. weekend flights from/to touristic places). On the other hand, airports lying in the 
intermediate range (20 ^ fc ^ 50) appear to follow a different distribution with stronger decay (a — 2). It might include airports 
which started or ceased to operate during the year and mainly airports which decreased the number of flights on a specific route 
or simply eliminated a route. 
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Figure 5: (Color online) Cumulative distribution of the number of flights divided by the number of routes {s\n+out / kin+out) for the BAN in 
2006. Data are separated on 3 intervals, each curve is normalized independently. The best fits are given by the functional form Cexp(— /3a;"): 
C = 0.803, /3 = 7.139x10"^ and a = 1.0 (for k ^ 10); C = 1.033, (3 = 5.409x10"^^ and a = 2.0 (for 20 s^ fc sC 50); and C = 2.004, 
f3 = 1.461x10"^ and a = 1.0 (for k ^ 70). 



Investigating the other two layers of the network, we observe that the traffic does not follow the shrink of the BAN. Indeed, 
the total traffic of passengers (layer 2) and cargo (layer 3) more than double in the period, from 18 millions passengers in 1995 
to 43 millions passengers in 2006, and from 0.4 megatons in 1995 to 0.9 megatons in 2006. Figure |6] shows that the average 
strength for each case (connections {s}nj^out)^ passengers {sin+out) ^^^ cargo (sf„^o„()) followed different dynamics, possibly 
converging to a higher occupancy level in the aircrafts (the capacity of the aircrafts increased as well). Important to observe 



that although this situation seems favorable to airline companies because of among other reasons, less maintenance costs (e.g. 
support staff or slots rent), actually, two large and old Brazilian airline companies ceased their activities 1 26] during this period 
and some of the routes were absorbed by other companies, two of them dominating most of the market nowadays. 
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Figure 6: (Color online) Average strength evolution on several layers in the BAN relative to year 2000: layer 1 representing the number of 
flights (s^); layer 2, the number of passengers (s^); and layer 3, the amount of cargo carried (s^). 

The analysis of the added and deleted connections in the BAN provides useful insights about the dynamics of flights (Fig.|7]i. 
For the sake of organization, we define eight categories: added connections being introduced between old airports on existing 
routes (aOOl)oron new routes (a002); added connections introduced between an old and a new airport (aON) and finally, added 
connections between two new airports (aNN), which is very unlikely (less than 0.01%). In the case of removal of connections we 
have the situation were they are deleted from old airports, keeping the route (dOOl) or removing the route all together (d002); 
the case where connections are deleted between an old and a removed airport (dOR), or deleted between two removed airports 
(dRR), which is as unlikely as case 4 (Fig.|7]i. Complementary to the routes analysis (Fig.|3]l, we observe that the BAN has large 
fluctuations in terms of the number of flights from year to year. Rather than trying to create new routes, companies follow the 
demand and invest on popular routes offering more flights. At the same time, connections tend to be removed before the route 
itself ceases to be operated. Moreover, it is reasonable to expect that when a new route is established, less frequent flights are 
offered to get a feedback before increasing the investment 12711 . 

50 



.if 



iHiliiill 



u 



_ 


aOOl 


nZi a002 


^ aON 


^ 


d001 


^ 


d002 


nn dOR 


^ 


aNN/dRF 



1 9961 9971 9981 99920002001 20022003200420052006 
Year 

Figure 7: (Color online) Percentage of added (a) (based in the current year) and deleted (d) (based in the previous year) connections between 
new (N), removed (R) and old (O) airports in the BAN. aOOl: added connections between old airports on existing routes; a002: added 
connections between old airports on new routes; aON: added connections between old and new airports; dOOl: deleted connections between 
old airports, but maintaining the route; d002: deleted connections between old airports, but removing the route; dOR: deleted connections 
between old and removed airports. aNN: added connections between new airports; and dRR: deleted connections between removed airports. 



These structural changes have strong impacts on single airports, with some of them becoming more central and others more 
peripheral with time. In terms of infra-structure, it means that the airports have different roles in the network and consequently. 





B 


f^in-\-out 




Airport 


1995 


2006 


A 


1995 


2006 


A 


1995 


2006 


A 


Manaus* 


0.110 (3) 


0.267 (1) 


+143% 


87 (11) 


85 (5) 


-2% 


50,222 (11) 


40,081 (11) 


-20% 


Sao Paulo* 


0.282 (1) 


0.254 (2) 


-10% 


205 (1) 


145 (1) 


-29% 


282,578 (1) 


346,495 (1) 


+23% 


Belem* 


0.070 (5) 


0.133 (3) 


+90% 


97 (8) 


72 (7/8) 


-26% 


62, 192 (6) 


47,995 (10) 


-23% 


Salvador 


0.041 (12) 


0.122 (4) 


+198% 


95 (9) 


95 (3) 


0% 


77, 822 (5) 


93, 757 (4) 


+20% 


Brasilia 


0.146 (3) 


0.085 (5) 


-42% 


164 (2) 


103 (2) 


-37% 


100, 449 (3) 


115,550 (3) 


+15% 


Belo Horizonte* 


0.061 (7) 


0.084 (6) 


+38% 


116 (4) 


81 (6) 


-30% 


61,925 (7) 


64,631 (7) 


+4% 


Rio de Janeiro* 


0.053 (8) 


0.077 (7) 


+45% 


118 (3) 


91 (4) 


-23% 


147, 165 (2) 


153,062 (2) 


+4% 


Curitiba* 


0.090 (4) 


0.035 (8) 


-61% 


111 (5) 


65 (10) 


-41% 


57, 674 (9) 


54, 507 (9) 


-5% 


Fortaleza 


0.069 (6) 


0.032 (9) 


-54% 


98 (7) 


72 (7/8) 


-27% 


61,473 (8) 


65,115 (6) 


+6% 


Recife 


0.012 (23) 


0.018 (14) 


+50% 


78 (14) 


63 (11) 


-19% 


81,813 (4) 


73, 386 (5) 


-10% 



Table II: The 10 most important cities in the years 1995 and 2006. The measurements correspond to betweenness centrality (B), total number 
of routes (kin+out) and total number of connections {sl„^o^t). In parentheses, the rank of the city in the respective year. "A" gives the 
variation based on the year 1995. Asterisks indicate cities served by 2 or 3 airports. 



policies related to new services and technologies have to be taken into account to avoid crashes of the whole system. The ten most 
important cities are selected and the respective centrality measurements are presented in Table HI] The number in parentheses 
shows the rank in the respective year and "A" gives the variation based on the year 1995. 

We note that some cities are served by more than one airport (indicated by asterisks), for instance, the city of Sao Paulo is 
mainly served by two airports ll28ll . Sao Paulo is consolidated as the major hub in the BAN, although suffering as one of the cities 
which lost more routes, it has the largest increase in the number of connections during the period. It suggests that some profitable 
routes motivate the increase in the number of connections, while economical uninteresting routes were removed. Although the 
betweenness decreased 10%, it is still one of the most central cities in the network ll29ll . Sao Paulo airports are important for 
several reasons: Sao Paulo is the main connection between Brazil and countries from other continents, which means that many 
people necessarily go there before traveling abroad (indeed, Sao Paulo is ranked as the 16th most central airport in the world 
airport network in the year 2000 |7]); the city itself is strategically located such that it connects the southern with the northern 
and the northeastern part of Brazil; and Sao Paulo is by far the most populous city in South America which demands a lot of 
traffic by all means of transportation. 

The interesting case of Manaus with a remarkable high betweenness is worth observing. This city is located in the middle of 
the Amazon rainforest (north of Brazil - see Figure [T] for the localization) and it has a central role to unite cities in that region 
with the other parts of the country. In general, most of the listed cities become more central with time (higher betweenness and 
more connections), indicating that though the network shrinks on the routes level, the increase in the amount of passengers and 
cargo is focused in specific profitable routes. 



IV. CONCLUSIONS 



Transportation networks have a strategic role in society because people and goods are constantly moving between distances 
of different scales. These movements are fundamental for diversity in culture and behavior, and also diffusion of information. 
Meanwhile, trades between different regions keep the economy strong, but require a high level of efficient, profitable and trustful 
transportation networks. On the other hand, transportation networks also contribute for instance, to disease spread, drugs traffic 
and smuggling. The aviation sector has a special relevance in this context since it provides fast, long-range paths between 
different places. 

Since the aviation sector is sensitive to economical fluctuations and geopolitical constraints, it is interesting to investigate how 
such network evolved during the recent years. By using the framework of complex networks theory, the structure and evolution 
of the BAN are investigated in several levels, from routes to passengers/cargo. The importance of the Brazilian airport network 
(BAN) lies in the fact that Brazil is an economically developing country with a large territory and sparse agglomerates of people. 

It is observed that although some structural features are in accordance with other previous results from different airport net- 
works, the BAN has a rewiring dynamics taking place at all levels. The results suggest that the companies have a tendency to 
invest in the most profitable routes rather than in new routes, consequently, increasing the number of connections on specific 
routes. The number of routes together with the number of airports decrease during the period, but the routes are constantly 
changing and not necessarily within the most connected airports. This dynamical evolution resulted on some airports becoming 
more central with time, while others become more peripheral. Specifically, the average betweenness of the network increased 
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during the period and the degree distribution become more homogeneous. Broad degree and betweenness centrahties (cumula- 
tive) distributions were observed in all years. Despite of the rewiring, the BAN showed a high number of cycles when compared 
to a random network, indicating that there are well defined internal structures on specific regions of the network. The structural 
changes in the ground levels had positive effects in the other levels, the number of passengers and the cargo carried through the 
network more than doubled in the considered period. 

Complementary questions have been left open and might be of interest. For instance, we do not know whether this structural 
evolutionary path is a pure effect of the socioeconomic development of the country or an universal feature of this type of system. 
Indeed, it seems that other countries/regions also have a dynamic network evolution, for instance, the number of available seats 
decreased in the US network and increased in Asia and Europe networks in the last 10 years ||23I1 . Moreover, specific system 
failures and attacks should be investigated more carefully since in fact, there is movement of carriers over these layers (e.g. 
congestion dynamics) and the respective time scale is different than the one investigated here (e.g. a delay occur within a period 
of one or few days). 
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[24] In this paper, new, old and removed airports correspond to new, old and removed nodes in the network. More specifically, old means that 
the node was in the network at least in the previous and in the current year. It does not directly correspond to constructed or demolished 
airports. 

[25] The code means: "a" or "d" for added (new) or deleted (removed) edges, and "N", "R" or "O" for new, removed or old airports respec- 
tively. 

[26] Indeed, one of the companies was sold and has undergone a restructuring process. 

[27] We might indicate possible external influences as the source of specific variations. For instance, the two most aggressive expansions in 
1996 and 1999 are apparently followed by large removal of connections in the next year, meaning that aggressive investments do not 
necessarily increases profits. The federal government changed the policy and removed market restrictions at the end of 1997 [181 . which 
clearly affected the sector as we can see in the expansion of the number of connections in the two following years (Fig.|7](. However, new 
regulation rules starting in 2003 |T8] can be the responsible for the large removal of connections in this year. The substantial removal 
of connections in cases 6 and 7 (in special) in the year 2005 can be related to problems with the previous cited airline companies which 
stopped operating. 

[28] Officially, Governador Andre Franco Montoro International Airport, which is in the neighbor city of Guarulhos and Congonhas Interna- 
tional Airport, which is located inside the city of Sao Paulo but no longer used as an international airport. 

[29] The evolution of Sao Paulo connections can be viewed in the support material. 



